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ABSTRACT 

We calculate Faraday rotation in global axisymmetric magnetohydrodynamic simulations of 
geometrically thick accretion flows. These calculations are motivated by the measured rotation 
measure (RM) of « —6 x 10^ rad m~^ from Sgr A* in the Galactic center, which appears 
to have been stable over the past ~ 7 years. In our numerical simulations, the quasi-steady 
state structure of the accretion flow, and the RM it produces, depends on the initial magnetic 
field threading the accreting material. In spite of this dependence, we can draw several robust 
conclusions about Faraday rotation produced by geometrically thick accretion disks: i) the time 
averaged RM does not depend that sensitively on the viewing angle through the accretion flow, 
but the stability of the RM can. Equatorial viewing angles show significant variability in RM 
(including sign reversals), while polar viewing angles are relatively stable if there is a large scale 
magnetic field threading the disk at large radii, ii) Most of the RM is produced at small radii for 
polar viewing angles while all radii contribute significantly near the midplane of the disk. Our 
simulations confirm previous analytic arguments that the accretion rate onto Sgr A* must satisfy 
Min ^ A/Bondi ^ 10~^ Mq yr~^ in order to not over-produce the measured RM. We argue that 
the steady RM w —6 x 10^ rad m~^ from Sgr A* has two plausible explanations: 1) it is produced 
at ~ 100 Schwarzschild radii, requires Min « 3 x 1O~*M0 yr~^, and we view the flow at an angle 
of ~ 30° relative to the rotation axis of the disk; in our simulations, the variation in RM across 
a finite-sized source is sufficient to depolarize the emission below w 100 GHz, consistent with 
observations. 2) Alternatively, the RM may be produced in the relatively spherical inflowing 
plasma near the circularization radius at ~ 10'^ — lO'' Schwarzschild radii, with the magnetic field 
perhaps amplified by the magnetothermal instability. Time variability studies of the RM can 
distinguish between these two possibilities. 

Subject headings: accretion, accretion disks - MHD - Galaxy: center 

1. Introduction cal modeling of stellar outflows and X-ray observa- 

tions of the diffuse thermal plasma in the central 
Observations of stellar orbits m the Galactic parsec show that winds from massive stars provide 

center indicate the presence of a « 4 x IO^Mq ^ mass supply of Meondi ~ 10"^ - lO-^M© yr"! 

black hole at the center of the Galaxy (e.g., ^t « 10^ Schwarzschild radii, which is approx- 

Schodel et al. || 2002t | Ghez et al. 11 20051 ). Numeri- imately the gravitational sphere of influence of 
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the black hole (iBaganoff et al. I I2OO3I: iQuataert 
20041 : ICuadra et al.l l2006h . The observed bolo- 



metric luminosity of the electromagnetic source 
Sgr A* coincident with the black hole is, how- 
ever, 5 orders of magnitude smaller than would 
be expected if all of this available gas were ac- 
creted with a radiative efficiency of 10%. As a 
result, the inflowing plasm a must have a very 
low radiative efhciency (e.g. jNaravan fc Yi Ill995r ) 
and/or very little of the gas supplied at large radii 
actually makes it down to t he black hole (e.g., 
iBlandford fc Begelman lll999f l. Global MHD sim- 
ulations of radiatively inefficient accretion flows 
show significant mass-loss consist ent wi th the lat- 
ter po ssibilitv Isto nc fc Pringle 2001 . hereafter 



SPOl; iHawlev. Balbus. fc Stone 1120011) . In addi- 
tion, local numerical simulations of electron heat- 
ing in coUisionlcss accretion flows flnd that the 
electrons receive a significant fraction of the grav- 
itational potential energy released in the accretion 
flow, requiring M « C AfBondi to explain th e low lu- 



minosity of Sgr A* (jSharma et al. 1120071) 



An important observational diagnostic of the 
physical conditions in the accreting plasma is po- 
larization measurements. Li near polarization from 



Sgr A* was first detected bv lAitken et al. 1 (|2000l) 



at sub-millimeter wavelengths; they found that the 
polarization fraction increases rapidly from 150 to 
400 GHz. Linear polarization is not detected at 
lower frequencies, although circ ular polarization is 



(e.g., iBower et al. Ill999l l2002t ). The decrease in 



the linear polarization fraction at low frequencies 
may be due to depolarization by di fferential Fara- 



day r otation across the source (e.g.. iBower et al. 
I999I ). In addition, the emission at low frequen- 



cies is optically thick, which suppresses the linear 
polarization if the par ticle distribution func tion is 
roughly thermal (e.g.. iGoldston etaL 2005 ). 

The detection of linearly polarized radiation 
at mm fc sub-mm wavelengths rules out high ac- 
cretion rate models in which M ^ -MBondi near 



the black hole (lAgol gOOOl : iQuataert fc Gruzinov 
2OOOI : iMeha et al.ll200Ci) . Analvtic estimates sho-w 



show 

that such models would produce a rotation mea- 
sure (RM) much larger than observed — so large, 
in fact, that all of the mm emission would be de- 
polarized by Faraday rotation. In the past five 
years, higher resolution constraints on the linear 
polarization of Sgr A* have generally confirmed 
the Aitken et al. results ( Bower et al. ] I2OO.4 



Marrone et al. 



Marrone et al. 



|2006| : iMacauart et al. I l2006f) . 
2007') performed the first mea- 
surement of the polarization of Sgr A* with suffi- 
cient sensitivity to determine the Faraday rotation 
in single-epoch observations, free from ambiguities 
introduced by possible time variation in the intrin- 
sic polarization angle. They found RM « — 6 x 10^ 
rad m~^ using observations at both 227 and 343 
GHz. The observed RM appears to be relatively 
constant over the 2 years of high sensitivity SMA 
observations (Marrone, private communication). 
Comparison with previous observations suggests 
that the RM towards Sgr A* has not changed sign 
over the past w 7 years. 

Propagation through the interstellar medium 
cannot account for the large RM observed towards 
Sgr A*. In addition. X-ray observations of the 
thermal plasma in the central parsec of the Galac- 
tic center shows that the electron number density 
is « 100 cm^'^ and the temperature is « 2 keV 



at a distance of « 1" 



0.04p c 



10^ 



fror 



the b lack hole ([Baganoff et al. 1 12OO3I: iQuataert 
20041 ). where Vg is the Schwarzschild radius. If the 



magnetic field on this scale were uniform, purely 
radial, and in equipartition with the thermal pres- 
sure, B K 2 niG and RM « 5000 rad m~^, two or- 
ders of magnitude smaller than what is observed. 
This shows that the observed RM must be pro- 
duced at distances ^ lO^r^ from the black hole, 
providing an important diagnostic of the physics 
of accretion onto Sgr A*. In view of the uniqueness 
of this probe, and the importance of Sgr A* as the 
paradigm for radiatively inefficient accretion flows, 
a more detailed study of the constraints imposed 
by the observed Faraday rotation is warranted. 

In this paper we carry out two-dimensional 
(axisymmctric) MHD simulations of accretion 
onto a central point mass and study Faraday 
rotation through the resulting turbulent accre- 
tion flow. This significantly improves on previ- 
ous analytic estimates of the Faraday rotation in 
radiatively inefficient accretion flows, which re- 
lied on simplifying assumptions such as power-law 
density profiles, uniform steady magnetic fields, 
etc. One caveat about our assumption of axisym- 
metry is that, by the anti-dynamo theorem, the 
magnetic field must eventually decay. However, 
the similarity of previous 2-D and 3-D simula- 
tions in the quasi-steady turbulent state, that 
lasts for hundreds of orbital periods at small radii 



2 



(jHawlev. Balbus. fc Stone 1 120011 ). gives us some 
confidence that our results will hold in 3-D. 

The remainder of this paper is organized as fol- 
lows. Our numerical techniques and methods for 
calculating RM are described in §2; we also briefly 
summarize the analytic scalings for RM in radia- 
tively inefficient accretion flows (i i2.3p . In §3 we 
discuss our results on Faraday rotation in MHD 
simulations. We focus on a fiducial model, but 
discuss the sensitivity of our results to the initial 
magnetic field configuration in the simulation, nu- 
merical resolution, and the size of the computa- 
tional domain. In §4 we summarize our results 
and discuss several scenarios that can account for 
the observed RM from Sgr A*. 

2. Methods and Analytic Scalings 
2.1. Numerical Simulations 

Unless specified otherwise, we use a numerical 
set up simflar to that of SPOl. We use the ZEUS- 

Tooiallbl ) 



2D MHD code (Stone k Norman 



with 

spherical (r, 0) coordinates to solve the equations 
of ideal MHD under the assumption of axisym- 
metry. The gravitatio nal potential is the pseudq - 
Newtonian potential of lPaczvnski fc Wiita l ( 1980l ). 
namely $ = -GM/{r - Tg), where Tg = 2GM/c^. 
The energy equation is adiabatic with the excep- 
tion of heating by shock-capturing artificial vis- 
cosity (which is a rather minor source of heating). 
We do not use an explicit resistivity to capture 
some of the energy lost to grid-scale averaging 
as it does not appear to significantly affect the 
dynamics of the flow (e.g., see SPOl). 

The initial condition is a geometrically thick, 
weakly magnetized, rotationally and pressure sup- 
ported constant angular momentum torus. The 
densit y in the initial torus is given by e quation 



(6) of IStone, P ringle. & BegelmanJ (199£) which 
is based on [Papaloizou & PringlCj (1984). The 
peak density of the torus is 1 in code units. The 
torus is truncated when the density falls below 
0.01. It is surrounded by a constant density corona 
(p 10-4); there is a density and temperature 
jump at the torus-corona boundary. 

The torus is specified by the radius of the ini- 
tial density maximum, tq, which we take to be 
ro = 40rg, lOOr^, or 200rg. The gravitational 
sphere of influence of Sgr A* on the surround- 
ing hot plasma is ^ lO^^g. Three-dimensional nu- 



merical simulations of accretion from stellar winds 
i mply a circulariza tion radius of ~ 10^ — lO^'rg 
( Cuadra et al ]l200()h . Ideally we would thus like to 
simulate an accretion disk which is ^ 10^ — lO^r^ 
in radial extent. This is, however, infeasible at 
the present time. Instead, our goal is to under- 
stand the generic properties of Faraday rotation 
through the steady state accretion flow structure 
that develops at radii < tq. 

Our grid is logarithmic with a higher density 
of points at small radii and in the midplane. 
The radial grid is spaced such that the num- 
ber of grid points from r,nin to (r,ninr„iax)^/^ is 
roughly the same as the number of grid points 

from (rmin^-max)^/^ tO Tmax, whcrC Tmin and Tmax 

are the inner and outer radial boundaries. Three 
different resolutions, 60 x 44, 120 x 88, and 240 x 
176, are used. We use outflow boundary condi- 
tions and set BrB^ < at rmin — 2rg. 

We consider two different initial magnetic field 
configurations in the torus: configuration A in 
which the field lies along the constant density con- 
tours (Fig. [ij and is derived from the vector po- 
tential — p^/Po with Po = 200 (the same as 
runs B fc F of SPOl), and configuration B which 
has a net radial field in the equatorial plane (see 
Fig. [8]& the Appendix). We find that the results 
of our Faraday rotation calculations remain sen- 
sitive to the initial field configuration even in the 
turbulent state. This is discussed in more detail 
in §3.2 and §4. 

The initial weakly magnetized torus is unsta- 
ble to the magnetorot ational instability (MRI; 
Balbus fc Hawlev]ll99ll ). which rapidly amplifies 
the magnetic field leading to MHD turbulence and 
efficient angular momentum transport. After an 
initial amplification stage, a quasi-steady turbu- 
lent accretion flow is established. However, at late 
times, due to the impossibility of a self sustained 
dynamo in 2D, the turbulent motions die away. 
We carry out the RM calculations in the quasi- 
steady turbulent state which lasts for hundreds of 
orbital periods at small radii. 

2.2. Rotation Measure 

Due to the different index of refraction for left 
and right circularly polarized radiation, the po- 
larization angle of linearly polarized light is ro- 
tated by an angle of A^RM on passing through 
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a non-relativistic magnetized plasma, where RM 
= {e^ /2TTmlc'^) J nB • dl is the rotation measure, 
e and nie are the electron charge and mass, c is the 
velocity of light in vacuum, n is the number den- 
sity of electrons, B is the magnetic field, and dl is 
the differential length along the line of sight. For 
a relativistic plasm a, RM is suppressed by a fac - 



tor of \n{9)/9'^ (jOuataert fc Gruzinov 
where 9 — kTe/uieC^. In addition, propagation 
through a relativistic plasma can lead to conver- 
sion between linear and circular polarization (e.g., 
Ruszkowski & Begelman 2002; Ballantyne, Ozel, 
& Psaltis 2007). We focus on non-relativistic plas- 
mas throughout this paper. 

We calculate RM along different lines of sight 
through the simulation domain in the quasi-steady 
turbulent state. The MHD variables {p, B, etc.) 
are output 40 times per orbital period at the den- 
sity maximum rg ; these variables are then used to 
calculate RM for different viewing angles. For Fig- 
ure [7] we used output 160 times per orbital period 
to guarantee that the short timescale variability 
in RM was captured. Our results are, however, 
largely insensitive to the exact time sampling used. 

The rotation measure in arbitrary code units is 
rescaled as follows 



code 



_ / (^^odc) 

0.002 



-3/2 



pB-dl (1) 



where (Mcodc) is the average mass accretion rate 
at Tniin — 2rg during the turbulent steady state, 
and Tin and rout are the inner and outer radii 
used in calculating the rotation measure. Numer- 
ical simulations with different resolution, tq, and 
initial field configurations have different values of 
{Mcode}- Equation ^ effectively re-scales all of 
the different simulations to have the same time 
averaged accretion rate, allowing a more useful 
comparison of the different simulations. The con- 
version of RMcodc to a physical rotation measure 
(in rad m~^) is then given by 



RM = 1.5x10 



10 



Mi, 



3/2 



BH 



10-5Moyr-i / \A x IO^Mq 

(2) 

where Mbh is the black hole mass and Mjn is the 
actual mass accretion rate through the horizon of 
the black hole (i.e., at rmin = 2rg). 



Physically rin in equation ^ corresponds to the 
radius where the electrons become non-relativistic, 
since there is little Faraday rotation from smaller 
radii. We choose Tin — lOvg unless specified other- 
wise. The true rjn in Sgr A* is probably somewhat 
larger, Tin ~ lOOr^ (e.g., Sharma et al. 2007), but 
we are restricted to smaller values of r^^ for nu- 
merical reasons. We discuss the scalings of our 
results to larger Hn in §4. The outer radius rout is 
chosen so that the density structure of the accre- 
tion flow is not affected by the initial conditions. 



We choose Tq 



30r„ for tq = 40r„, rout = 60rg 



for ro = 100r„, and ro 



120r„ for ro = 200r„ 



We calculate the RM in two ways in the sim- 
ulations: along individual lines of sight {9 = con- 
stant in spherical coordinates) and in a cylindri- 
cal "beam" of diameter I?bcam- Physically, the 
relevant beam size is expected to be the size of 
the synchrotron source which is itself a function 
of frequenc y, with Ubeam ^ 10ro(i//100 GHz) ~^ 



cy . , , 

for Sgr A* (iBower et al. ll2004H Shen et al. II2005I: 



Loeb fc WaxmanI |2007| ) We take D, 



beam 



lOrg 

but find that our results are not sensitive to a fac- 
tor of ~ 2 variation in this choice. 

For an individual line of sight we evaluate the 
integral in equation ([T]) by adding up the contri- 
bution of each 9 = constant grid point that lies 
between rjn and rout • For a beam of width -Dbcam 
all grid points lying within the beam contribute 
to the RM integral. This effectively assumes that 
the source has a constant surface brightness across 
the synchrotron photosphere. Since the size of the 
grid cells is not uniform, we calculate the beam 
averaged RM by calculating the average value of 
nJ5|| within the beam, weighing the contribution 
of each grid cell by its area, and multiplying by 
''out — ?'in (the radial length of the region contribut- 
ing to the calculated RM) . Using the RM for each 
time-slice in the turbulent state we calculate the 
time averaged RM and the standard deviation in 
time, for both individual lines of sight and beams. 
The standard deviation in time provides a useful 
measure of the time variability of the RM. We also 
calculate the spatial variation in RM within the 
be£fftf'fa| width I?boam; this can be used to quan- 
tify the' effects of beam depolarization. 

2.3. Analytic Scalings 

To provide some context for the numerical sim- 
ulations that follow, we briefly review the an- 
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alytic expectations for RM in radiatively ineffi- 
cient accretion flows. For a flow with a density 
profile of n oc and a virial tempera- 

ture T oc , the equipartition magnetic field 
strength should vary as i? oc {nTY^'^ oc r~^/**+P/^. 
Thus the local contribution of a given decade 
in radius to the rotation measure is given by 
dRM/dh-ir oc rnB oc r-'^/4+3p/2. For Bondi ac- 
cretion, p — and dRM/dhir oc r^^/*; RM is 
thus dominated by small radii. For Sgr A*, the 
Bondi scalings may be applicable o utside the cir- 
cular ization radius of ^ 10'^ — lO'^rg (ICuadra et al 
20061 ). although this need not be true in the 



presence of magnetic fie lds or significant heat 
conduction (e.g., J gumcns hchev fc Naravan 20021 : 
Johnson fc Quataertll2007l ). 

In contrast to the B ondi model, analytic models 

(iBland ford fc Begelman 1999:,Quataert fc Gruzinov 



2001 



2000 ) and numerical simulations (IStone fc Pringle 



iHawlev. Balbus. fc Stone 2001 ) of rotating 
radiatively inefficient accretion flows find p 1, 
which corresponds to an accretion rate which 
decreases with radius as M oc r. In this case, 
dRM/d\nr oc r~^^'^. There is thus a significant 
contribution to RM from both small and large 
radii in the flow. 

3. Numerical Results 

3.1. The Fiducial Simulation 

As a fiducial case, we focus on the simulation 
with tq = 200rg, a resolution of 120 x 88 grid 
points, and with initial magnetic field unit vec- 
tors shown in Figure [1] ( "configuration A" ) . The 
resolution in the innermost portions of the torus 
is such that the fastest growing MRI mode is re- 
solved with about 5 grid points. Different resolu- 
tions, ro, and initial magnetic field configurations 
are discussed later. Several properties for these 
simulations are summarized in Tables [T] fc [2l 

The shear in the initial torus converts radial 
magnetic field into azimuthal field. Simultane- 
ously, the MRI amplifies the magnetic field. The 
dynamical timescales are shortest at small radii in 
the torus, where mass flows in as it loses angular 
momentum because of magnetic stresses. Figure [5] 
shows the mass accretion rate through = 2rg 
as a function of time for this simulation. A quasi- 
steady turbulent accretion flow is established after 
5 orbital periods at 100 Vg (note that 1 orbital pe- 



riod at 100 rg is « 3 (AfBH/4 x lO'^Mg) days). The 
accretion rate is, however, highly variable and can 
sometimes change by an order of magnitude on 
short time scales (of order the orbital timescales 
at small radii). This is somewhat reminiscent of 



the fl aring activity seen in Sgr A* (jBaganoff et al. 
20011 : iGenzel et al. Il20n3h . 



Figure [3] shows contour plots of the density and 
plasma /3 (the ratio of thermal to magnetic pres- 
sure) for the fiducial simulation, averaged over the 
quasi-steady turbulent state. Poloidal magnetic 
field unit vectors are superimposed on the den- 
sity contours. Figure [3] shows that a high density, 
moderately thin disk with very high f3 is formed in 
the equatorial plane of the initial torus (similar to 
Figures 4 fc 5 of SPOl). The accretion flow in the 
turbulent state is influenced by the initial mag- 
netic field configuration. Shearing of the initially 
radial magnetic field above and below the equa- 
tor leads to a strong toroidal field in the corona 
which compresses the midplane into a dense mod- 
erately thin disk with high /3. By contrast, both 
the density and P are quite small in the polar re- 
gions. The average magnetic field vectors shown 
in Figure [3] for the quasi-steady structure clearly 
maintain memory of the initial magnetic field con- 
figuration shown in Figure [T] 

For our fiducial simulation, we use rjn = lOj'g 
and Tout = 120rg as the inner and outer radii in 
the RM integral (eq. [T]). Figure 2] shows the 
time averaged RMcode and the standard deviation 
in time along different lines of sight. As expected, 
the RM is fairly symmetric about the midplane 
[9 — 90°), except for the sign change which arises 
because of the initial field configuration. In con- 
trast to (RMcodc), the standard deviation peaks 
at the midplane and is significantly smaller near 
the pole. This is because the midplane is highly 
turbulent, while the polar regions contain a com- 
paratively stable magnetic field. Figure |4] shows 
that the time variability in RM is larger than the 
average RM except within « 40° of the poles. 

As discussed in ^2.2[ a more apt comparison 
to observations is the Faraday rotation averaged 
over a finite source size, rather than along in- 
dividual lines of sight. Figure O shows RMcode 
and the standard deviation as a function of 9 for 
-Dbeam = lOfg. The results in Figure[5]are reason- 
ably similar to the single angle results shown in 
Figured! The primary difference is that the fluc- 
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tuations in the RM in the midplane {9 ^ 90°) are 
smaller by a factor of ~ 3 because the beam aver- 
ages over the spatial fluctuations in the turbulent 
disk. Nonetheless, the temporal fluctuations im- 
plied by Figure [5] are still sufficiently large that one 
would only expect a relatively steady RM for polar 
viewing angles. Note also that both {RMcode) and 
its standard deviation do not depend that strongly 
on viewing angle 0; there is roughly a factor of ~ 5 
variation from pole to equator. 

To quantify which radii dominate the contribu- 
tion to the RM, it is useful to consider the log- 
arithmic derivative dRM / d\nr. Figure [6] shows 
dRMcodc/d In r as a function of radius for an equa- 
torial and 6 = 30° beam with i^bcam = lO^gS 
both the time average and standard deviation are 
shown. dRMcodc/dlnr fluctuates with radius in 
Figure [6] because we are taking a numerical deriva- 
tive with a finite number of grid points in the 
beam. For the equatorial line of sight, both the 
time average (solid) and standard deviation (dot- 
dashed) of RM receive roughly equal contributions 
from all radii. This is consistent with the analytic 
scalings reviewed in §2.31 In addition, as already 
shown in Figure^ the standard deviation is signif- 
icantly larger than the time averaged RM for the 
equatorial beam, i.e., the RM is highly variable in 
the turbulent disk. By contrast, the 6 = 30° beam 
has a (dRMcodc/d^nr) (dashed) which is some- 
what larger than the standard deviation (dotted) . 

For the polar line of sight through the corona 
of the disk, both {dRM code / dhir) and the stan- 
dard deviation are dominated by small radii, in 
contrast to the equatorial lines of sight where all 
radii contribute significantly. The solid line super- 
imposed on the numerical results in Figure [S] is a 
dRM code/ din r ^ r^-^/"^ power law which does a 
reasonable job at fitting the numerical results. 

Spatial variations in the amount of Faraday 
rotation across a finite-sized source can lead to 
depolarization of the source at wavelengths where 
a±{RM)\'^ > 1, where aj_{RM) is the spa- 
tial standard deviation in R M across the source 
beam depolarization" ; e.g., iBurn lligedlTribble 



199ll ). To quantify this effect in our simula- 
tions. Figure [7] shows the beam averaged RM and 
a±{RM) as a function of time for 6 = 30° and 
9 = 90°, taking Z^beam = lOrg and ^ lOr^; 
a±{RM) is calculated by subdividing the lOr^ 
wide beam into smaller beams. The time aver- 



age value of a±{RM)/RM is « 10 for 9 = 90° 
and « 0.3 for 9 = 30°; it is somewhat smaller for 
9 < 30°. For a larger value of rin = 25rg we find 
similar results for 9 — 90° but find that the spatial 
variation in RM is somewhat larger at 6* = 30°, 
with an average value of a±{RM)/ RM « 0.9. 

3.2. Effects of the Initial Magnetic Field 
Configuration 

The results of the previous section show that 
several properties of the quasi-steady turbulent 
state in the accretion flow at late times are deter- 
mined by the initial magnetic field threading the 
torus at the beginning of the simulation. In par- 
ticular, the antisymmetry of RM about the equa- 
torial plane (Fig. [5]) and the relatively coherent 
magnetic field at polar angles (Fig. [3] and [5]) are 
all signatures of the initial magnetic field struc- 
ture shown in Figure [TJ It is thus important to 
understand which of our results on the disk struc- 
ture and the RM are robust to changes in the ini- 
tial field configuration. Figure [8] shows the initial 
magnetic field unit vectors for the alternative ini- 
tial field configuration we considered ("configura- 
tion B" ; see the Appendix for the vector potential 
used to generate this field). In this case, the initial 
field is symmetric about the midplane and has a 
net radial magnetic field in the equatorial plane 
of the initial torus; /3 is initially w 10 — 30 in the 
equatorial plane and « 100 — 10^ in the rest of 
the torus. Note that although the initial magnetic 
field is still confined to the torus, the field lines are 
not along the constant density contours (as they 
were in Fig. [1]). In addition, the vertical scale 
of the magnetic field structure in configuration B 
is roughly half that of configuration A. Thus we 
need a higher resolution to resolve the dynamics of 
the MRI. We take the simulation with rg = lOOr^ 
and a resolution of 240 x 176 as the standard case 
for initial field configuration B. For this case the 
fastest growing MRI modes at 120rg (just above 
and below the equator) are resolved with about 8 
grid points. 

Figure [9] shows 2-D contour plots of the density 
and plasma /3 averaged over the turbulent steady 
state. The resulting disk is significantly thicker 
compared to that arising from field configuration 
A (see Fig. [3]). In addition, the disk scale-height 
is roughly proportional to the radius as in analytic 
models of radiatively inefficient accretion flows. 
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Unlike configuration A where l3 ^ 1 in the disk 
midplane, here P ^ 1 in the turbulent disk. The 
average field structure in the quasi-steady state is 
again reminiscent of the initial field configuration 
with radially inward field lines in the equatorial 
plane. The thick disk is surrounded by a rela- 
tively high /3 turbulent region, with fairly low /3 
plasma at the poles. In general, the density and 
magnetic field structure of the disk is rather dif- 
ferent from that in Figure O This is also shown 
in Figures [10] & [Til which compare the density 
and (3 within 12° of the equatorial plane averaged 
over the turbulent state as a function of radius for 
the two different initial magnetic configurations. 
The magnetic field is stronger for configuration 
B, and f3 decreases with radius, likely because of 
flux freezing of the initial radial field. In addi- 
tion, the density does not increase significantly 
inwards for configuration B. These differences in 
the steady state disk structure appear be a direct 
consequence of the different initial magnetic field 
structure. Although a net inward radial field may 
not always be a realistic initial condition, part of 
the motivation for choosing this particular config- 
uration was to understand whether the RM in the 
bulk of the disk showed significant time variabil- 
ity even in the extreme case of a coherent initial 
radial field. 

Figure [12] shows the time averaged RM and the 
standard deviation with I?bcam = 10»'g as a func- 
tion of angle 6. The RM is roughly symmetric 
about the equatorial plane and peaks at 6* « 90°. 
This is a direct consequence of the initial magnetic 
field configuration, which is symmetric about the 
equator and has a net radial field in the equator of 
the torus. Figure [T^ also shows that the time vari- 
ations in RM are larger than (RMcodc) at essen- 
tially all angles. This is in contrast to the results 
for initial magnetic field configuration A shown in 
Figure [5] where (RMcodc) is larger than the fluc- 
tuations near pole. In the present case, the polar 
viewing angles are relatively turbulent and show 
magnetic field reversals in time. Interestingly, al- 
though (RMcodc) is non-zero in the equator, the 
fluctuations due to the turbulence are still suffi- 
ciently large that one would expect signiflcant time 
variability and sign changes in RM for equatorial 
viewing angles. 

The time averaged dRMcodc/dlnr and its stan- 
dard deviation for equatorial and 6 = 30° beams 



with Dboam = lO^g are shown as a function 
of radius in Figure 1131 For the polar viewing 
angle, the smallest radii dominate the contribu- 
tion to both (RMcode) and its time variation, 
with dRMcodc/dlnr cx being a reasonable fit. 
For the equatorial beam, however, both (RAIcodc) 
and the standard deviation are dominated by the 
largest radii. This is in contrast to the results 
shown in Figure [6] for initial field configuration A, 
and in contrast to the simple analytic estimates 
from The significant contribution of large 

radii to (RMcodc) in the equatorial plane is in 
part because the density in the equator increases 
with radius for initial field configuration B (see 
Fis.lT 



3.3. Convergence Studies 

To test whether our results depend on the res- 
olution of the simulation, we carried out simula- 
tions at a number of different resolutions: 60 x 44, 
120 X 88, and 240 x 176 for initial configuration 
A (see Table P, and 120 x 88 and 240 x 176 for 
initial field configuration B (see Table[21) . The den- 
sity and magnetic field structure in the turbulent 
state is similar for all resolutions with initial con- 
figuration A. However, configuration B requires a 
higher resolution to resolve the dynamics of the 
MRI. As a result, the lower resolution simulation 
(120 X 88) has a smaller (M) than the higher res- 
olution simulation (240 x 176). 

As a dimensionless measure of whether our re- 
sults for the RM have converged, we define the 
coherence parameter (Coh) as 



Coh = 



-'^in 



nB ■ dl 



nlB -dll 



(3) 



Coh quantifies the effect of cancellation due to 
magnetic field reversals on the inferred RM. If the 
large-scale dynamics of the magnetic field is prop- 
erly resolved, one would expect that simulations 
with different resolution would have statistically 
similar values for Coh. 

Figure [14] shows the beam averaged Coh at 
9 = 90° (the equator) as a function of time for 
initial configuration A for three different resolu- 
tions and ro — lOOr^. All three resolutions show 
statistically similar variations about (RM) k, 
(corresponding to a similar average and standard 
deviation for RM; see Table [1]). The time variabil- 
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ity is also quite similar, implying that the resolu- 
tion does not appear to be significantly affecting 
our determination of RM. 

Figure shows the beam averaged Coh for ini- 
tial configuration B as a function of time in the 
equatorial plane for tq — lOOrg. The results are 
similar for the two different resolutions (240 x 176 
and 120 x 88). RM is negative at most times be- 
cause of the initial magnetic field. There are, how- 
ever, times when RM is positive, consistent with 
the large standard deviation shown in Figure [121 

Ideally for Sgr A*, one would like to simulate 
an initial torus at tq 
ner disk extends to ~ 
expected circulariz ation radius in Sgr A* (e.g., 
Cuadra et al. l2006h . Because of finite computa- 
tional resources, however, we can only simulate 
a smaller disk. We carried out simulations with 



~ lO^rg, so that the in- 
10^ - lO^'rg, similar to the 



ro = 40rg (rout = 30rg), lOOrg (ro 



60rg), & 



200ro 



120rg) for both initial field config- 



urations A & B (see Tables [T] & [2]). Although 
there are some quantitative differences between 
the different simulations, the basic results for RM, 
its statistical fluctuations, and its variation with 
radius, do not seem to depend significantly on 
the radial extent vq, particularly when compar- 
ing the To = lOOrg and ro = 200rg simulations 
(the ro = 40rg simulation has comparatively little 
radial dynamic range). This provides some confi- 
dence in the numerically calculated properties of 
RM through the accretion fiow. 

4. Discussion 

We have carried out an analysis of Faraday ro- 
tation in global axisymmetric MHD simulations 
of radiatively inefficient accretion flows, motivated 
by the observations of linear polarization and 
Faraday rotation in Sgr A*. We first summarize 
our primary results and then discuss their impli- 
cations for models of Sgr A*. 

We find that the accretion disk structure, and 
in particular the rotation measure RM, retains 
memory of the initial magnetic field threading the 
plasma, even in the quasi-steady turbulent state at 
late times. Physically, this is equivalent to a de- 
pendence on the properties of the magnetic field 
in the plasma feeding the accretion flow at large 
radii (which is typically hard to determine obser- 
vationally). We considered two very different ini- 



tial magnetic fleld geometries, one with a net ra- 
dial field in the equator of the initial torus and one 
without (see Figs. [S] [configuration B] &[T] [config- 
uration A], respectively). In configuration B, the 
initial magnetic field unit vectors are symmetric 
about the equatorial plane, while in configuration 
A they are antisymmetric. The initial symmetry 
of the magnetic field remains in the RM as a func- 
tion of viewing angle 9 determined in the turbulent 
accretion flow (compare Figs [5] & [T^ . In addi- 
tion, the steady state density and f3 proflles are 
rather different for the two different initial mag- 
netic field conflgurations (Figs. [TUlfclHI). Given 
the anti-dynamo theorem, it is not that surpris- 
ing that the initial field structure imprints itself 
on the late-time disk structure. However, the sim- 
ilarity of previous 2D and 3D global d isk simu- 
lations (jHawlev. Balbus. &: Stone 1120011 ) suggests 
that this will also be true in three dimensional 
simulations. It remains to be seen if simulations 
of accretion disks in three dimensions lose memory 
of the initial conditions at very late times (many 
orbital periods at large radii), e.g., by undergoing 
an inverse cascade which generates a large-scale 
magnetic field, or if the dependence on the large- 
scale initial magnetic field that we see is a generic 
and robust feature of accretion disk dynamos. 

In spite of the dependence on the initial mag- 
netic field threading the plasma, we can draw a 
number of general conclusions about the proper- 
ties of Faraday rotation though magnetized geo- 
metrically thick accretion fiows. 

1. The RM through the midplane of the tur- 
bulent disk is highly time dependent, with a 
typical dispersion in RM that is larger than 
the mean. This implies significant changes 
in the observed sign of the RM. We see sign 
changes in RM even for our initial condition 
with a net radial field in the midplane of the 
torus (configuration B; Fig. [T^ . an initial 
condition that was chosen to maximize the 
likelihood of a coherent sign of RM. 

2. In the turbulent disk, RM receives roughly 
equal contributions from all radii (Fig. [S] 
& [T5|) . Analytic scaling arguments for disk 
models that are consistent with the simula- 
tions also suggest a weak dependence of RM 
on radius, with small radii dominating by a 
modest amount {dRAI / d\nr oc r~^/*; 
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3. RM is only a modest function of viewing 
angle 9 (see Figs [5] & [E]); the polar RM 
is typically ~ 0.2 of the equatorial RM for 



nn = lOrg. 



4. If there is a large scale magnetic field thread- 
ing the plasma at large radii (our configu- 
ration A), the resulting corona of the disk 
also has a larg e-scale, stable ma gnetic field 
(see also SPOl.I Pe Villiers II2006I ). The RM 
viewed through the corona has a dispersion 
in time smaller than the mean (Fig. O and 
is thus relatively stable. 

5. For polar viewing angles, the RM is primar- 
ily produced at small radii; most of the con- 
tribution will thus arise from near the ra- 
dius at which the electrons become relativis- 
tic (Fig. El & [131). 

6. There is significant spatial variation in the 
amount of Faraday rotation through a finite 
sized beam, with (j±{RM) ~ 0.3 RM at po- 
lar viewing angles (Fig. [7]). This spatial 
variation in RM can cause depolarization at 
frequencies for which (t±{RM)X'^ > 1. 

4.1. Implications for Sgr A* 

Our resul ts con f irm previous analytic argu- 



-P 

ment s (e.g., lAgol I l2000t iQuataert fc Gruzinov 



2000l ) that models of Sgr A* with M - Msondi ~ 
10~^ Mq yr~^ are ruled out by the measured 
Faraday rotation of RM sa —6 x 10^ radm^^. For 
an accretion rate of M ~ AfBondi, the electrons 
must be marginally relativistic even at small radii 
in order to account for the low luminosity from 
Sgr A*. In this case, our choice of Hn = lOvg 
as the minimum radius for which Faraday rota- 
tion is produced (see eq. [Tj) is quite reason- 
able ( and even conservativ e; in the ADAF mod- 
els of iNaravan et al.lll998l the electron tempera- 
ture never exceeds 6 x 10^ K). Our simulations 
show that for nearly all viewing angles 9 and 
for both of the initial magnetic field configura- 
tions we considered, the dimensionless rotation 
measure RM^odc > 10"^ (Figs. E]&[121). Even 
though the time averaged RM may be smaller 
than this for certain viewing angles (e.g., through 
the equator of the turbulent disk, 9 — 90°, in 
Fig. [5|), the time variations are significant and so 
RMcodc 10~^ is a good indication of the typical 



instantaneous RM that would be observed at any 
9. Using equation ^ to convert RM to real units, 
we find that M Afeondi ~ 10~^ Mq yr~^ would 
imply RM ^ lO^*' radm^^, over 4 orders of mag- 
nitude larger than the observed value. Thus such 
models are strongly ruled out. 

We now discuss two scenarios that can account 
for the observed RM « —6 x 10^ radm~^ and its 
apparent stability over the past ~ 7 years. These 
scenarios each identify the RM-producing region 
with a key radial scale in the inflowing plasma. 

First, our simulations show that if there is a 
relatively coherent magnetic field in the plasma 
at large radii (configuration A) and if the ob- 
server's line of sight is through the polar regions 
of the accretion flow {9 ^ 30°), then the time 
variability of RM appears to be sufficiently small 
to account for stable sign of RM observed from 
Sgr A* (Figs. E] & E]). For these polar viewing 
angles, we find RMcodc ~ 3 x 10~^, i.e., RM « 
4.5 X lO^" (Min/10-5 Mq yr-i)3/2 radm"^. The 
RM at polar viewing angles is produced at small 
radii with dRM/dlnr oc r~" and awl — 1.5. 
Our simulations have = lOr^, but lOOrg 
is probably a more accurate estimate of the radius 
at which the electrons become no n-relativistic in 
Sgr A* fe.g.. lSharma et al.1l2007l ). Rescaling our 
simulations to a larger value of rin yields RM w 
10^ (Min/10-5 Mq yr-i)3/2 (ri„/100rg)-i-5 radm'^ 
for a = 1.5. Thus RM « 6 x 10^ radm^^ 
requires Mn « 7 x 10"* (nn/lOOrg) Mq yr'^. 
For a w 1, which is also reasonably consis- 
tent with our simulations, we find Mjn « 3 x 
10~* (nn/100rg)2/3 M© yr-^ Because ri„ is only 
moderately well-constrained, this estimate of Mjn 
is probably uncertain by a factor of ~ 5. In spite 
of this uncertainty, our estimated Mjn required 
to account for the observed RM from Sgr A* is 
comparable to previo us analytic estimates (e.g., 



Marrone et al. I l2007f ): it is also similar to the 



accretion rate estimate of [Sh arma et al. I (|2007l ). 
who inferred Min by calculating the radiative effi- 
ciency of the accreting plasma in Sgr A*. 

A second possible explanation for the observed 
RM from Sgr A* is that it is not directly produced 
within the nearly Keplerian accretion flow at small 
radii, but is instead produced in the more slowly 
rotating plasma further from the black hole. Nu- 
merical simulations of accretion onto Sgr A* from 



stellar winds located at ^ 10 



lO^rg find that the 
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circul arization radius of t he flow is Tcirc ^ 10^ — 
lO^rg (jCuadra et al.ll2006l ). As discussed in §1, the 



RM produced by the observed thermal plasma at 
~ lO^rg is « 5000/3"^/^ radm~^, where /3 is the 
ratio of the thermal energy density to the magnetic 
energy density, assuming a magnetic field with a 
significant radial component. If the Bondi solu- 
tion is applicable outside rcirc, then dRM/dlnr oc 
r~^/"' and the RM produced by the roughly spher- 
ical inflowing plasma at large radii is RM 



10^/3 



-3/2 I 



c/lO'^rg) radm ^ (roughly in- 



dependent of viewing angle). Circularization radii 
in the range expected can thus plausibly account 
for the observed RM from Sgr A* if /3 ^ 1. 
Note that in this interpretation, Mjn must be 
< 7 X 10"® (nn/lOOrg) Mq yr-i (for a = 1.5) 
in order for the Faraday rotation from r <g; rcirc to 
be small compared to that from ~ ^circ- 

This interpretation requires a significant mag- 
netic field at large radii, with ~ 1. Such a field 
could perhaps be produced by flux freezing of a 
moderately weaker field residing in the central 
parsec ~ lO^r^. More interestingly, spherically 
inflowing plasma in Sgr A* is unstable to the 
magnetothermal instability (MTI) which ampli- 
fies magnetic fields in a low collisionality plasma 
when the temperature increases in the directio n 
of gravity (|Balbus Il200ot IParrish fc Stone Il2005h . 
The MTI preferentially operates at large radii in 
spherical accretion flows (> lO'^rg) because only 
there is the growth time sh ort compared to the 
inflow time (e.g.. Fig. 9 of Johnson fc QuataertI 
[2007). In a hydrostatic atmosphere with a flxed 
temperature difference between the top and bot- 
tom of the atmosphere, the MTI leads to a signifl- 
cant radial heat flux, magnetic field amplification, 
and field lines combed o ut i n the radial direction 
Parrish fc Stone 2005 ). Johnson fc QuataerO 



20071 ) argued that the heat fiux associated with 



the MTI could modify the dynamics of spherical 
accretion at large radii, which would in turn alter 
dRM/dlnr from the Bondi estimate above. Thus 
the efficacy of the MTI for producing the observed 
Faraday rotation from Sgr A* needs to be evalu- 
ated in a self-consistent dynamical calculation. 

The above interpretations of the Faraday ro- 
tation observed towards Sgr A*, namely that it 
arises either from rin ~ lOOr^ or rdrc ~ lO^'^gi can 
probably be distinguished by careful monitoring of 
the variability of RM, which should be present at 



some level even in the absence of sign changes. In 
the Tin scenario one would expect variability on ~ 
hour-day timescales, while in the rdrc scenario, the 
variability would occur on much longer timescales, 
of order years. The apparent lack of sign i ficant 



RM v ariability to da te (e.g.. Bower et al. 2005t 



arrone et al. Il2007| ) may favor models in which 
the RM is produced at large radu ^ rdrc- 

If radii ~ rjn contribute significantly to the 
observed RM, our simulations can also be used 
to quantify whether Faraday depolarization can 
explain the rapid drop in the polarization frac- 
tion observ ed at frequencies b e low a few 100 



GHz (e.g., iBower et al. I 120031 : iMarrone et al. 
20061 ). Beam depolarization becomes important 



when (7±{RM)X^ > 1, where a± is the vari- 
ation in RM across a finite-sized source (e.g.. 
Burn 1 ll966HTribbleiri99F ). Our simulations have 
a±{RM) ~ 0.1 — 1 RM for polar viewing an- 
gles, with the exact value depending on 9 and rjn 
(Fig. [7] and gST]). Taking aj_{RM) « 0.3 i?A/ 
as a fiducial value, and using the observed RM of 
« 6x lO^rad m~^, this implies that emission below 
« 150 GHz should be depolarized by beam depo- 
larization, in reasonable agreement with observa- 
tions. Because of the variations in a± in our cur- 
rent simulations with r-^ and 6, higher resolution 
simulations and simulations with rjn ^ lOOr^ will 
be required to more accurately determine the fre- 
quency below which beam depolarization becomes 
important. Figure [7| also shows that a±{RM) in- 
creases by a factor of ~ 10 at certain times. This 
could give rise to a significant decrease in the ob- 
served polariza t ion si milar to that observed by 



Marrone et all (|2006l ) at 340 GHz. 



In a recent paper, iBallantvne et al.l ([20071) ar- 
gue that the "normal" Faraday rotation we have 
focused on in this paper may not be the dominant 
process shaping the observed polarization from Sgr 
A*. Instead, they argue that because the plasma 
is relativistic over a large range of radii, the nor- 
mal modes of the plasma are elliptical. In this 
case, the observed polarization probes the proper- 
ties of the normal modes in the radio-emitting re- 
gion, which in turn depends on the magnetic field 
and temperature structure of the plasma. How- 
ever, current observations of Sgr A* show that the 
change in the polarization angle from 100-400 GHz 
is consistent with the variat ion expected for 
normal Faraday rotation (e.g., iMacquart et al. 
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20061: iMarrone et aTl l2007t ). This strongly sug- 



gests that non-relativistic plasma at large radii is 
modifying the intrinsic polarization of the source 
via Faraday rotation. It also requires that the 
intrinsic polarization angle of the source is rela- 
tively independent of frequency. Finally, we note 
that the source size is inferred to decrease signifi- 
cantly at the frequencies for which linear polariza- 
tion is detected, from ss lO rn at 1 00 GHz to « 2rg 



at 400 GHz ( Bower et al. 1120041 ). Thus different 



frequencies sample different plasma along the line 
of sight propagating to observers on Earth. To 
produce the same RM independent of frequency 
requires that the source size at all frequencies is 
much smaller than the size of the region where 
the RM is produced (so that there is no variation 
in RM with frequency). This criterion is reason- 
ably satisfied for both of the interpretations of the 
Faraday rotation from Sgr A* considered above. 

Throughout this paper wc have focused on the 
implications of our results for Sgr A* because of 
the wealth of observational data available. We an- 
ticipate, however, that as observational techniques 
improve our simulations will prove useful for in- 
terpreting the physical conditions in the vicinity 
of other accreting black holes. 
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A. Vector potential for configuration B 

To define the vector potential for initial field configuration B, we set up a new cylindrical coordinate 
system whose origin is on the equator of our standard spherical coordinate system at a distance of 1.3ro 
(roughly the mid-point of the initial torus), and whose x and y axes are aligned with the local —6 and —r 
axes of the spherical coordinate system, respectively. The vector potential for initial field configuration B is 
then given by 

for Tcy < Tm, 

otherwise, (Al) 

where fey is the radius in the cylindrical coordinate system, = 0.5ro, and is the angle with respect to 
the X- axis in the cylindrical coordinate system. The vector potential is located at cell corners and can be 
differenced to get the magnetic field unit vectors shown in Figure [H 



A = 0.005 cos ( ^-^ ) cos((?!)) 
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Table 1 

RM FOR Initial Field Configuration A 





Resolution 


(Mcodc) 


(Coh)" 


<T(Coh) 






(Coh) 


(T(Coh) ^ 




0-(i?M,„de) 


200 rg" 


120 X 


88 


0.0025 


-0.1 


0.22 


-1.1 xlO~* 


4.1 xlO~* 


0.85 


0.13 


2.8 xlO~* 


10-* 


200 rg 


60 X 


44 


0.002 


0.04 


0.16 


8.3 xlO"-"* 


2.7 xlO"" 


0.95 


0.088 


2.1 xlO~* 


10-* 


100 T-g 


120 X 


88 


0.0029 


-0.05 


0.14 


-1.6 xlO"'' 


3.9 xlO"'' 


0.83 


0.14 


1.6 xlO~* 


5.7 xlO-^ 


100 r. 


60 X 


44 


0.0022 


-0.04 


0.17 


-8.5 xlO"'^ 


2.6 xlO"'' 


0.93 


0.12 


2.3 xlO""" 


6.8 xlO"^ 


100 rg 


240 X 


176 


0.0033 


0.04 


0.18 


8 Xl0~^ 


3.4 xlO"* 


0.63 


0.26 


1.4 xlO~* 


9 xlO-^ 


40 rg 


120 X 


88 


0.0031 


-0.06 


0.12 


-2.2 xlO~* 


3.7 Xl0~* 


0.66 


0.24 


9.5 xlQ-'^ 


5.7 xlQ-^ 


40 rg 


60 X 


44 


0.0028 


-0.21 


0.17 


-4.2 xlO-* 


3.1 XIO"* 


0.84 


0.15 


1.9 xlO-" 


8.1 xlO-^ 


40 rg 


240 X 


176 


0.0039 


-0.03 


0.1 


-1.4 xlO-"" 


4 xlO-"" 


0.69 


0.22 


8.5 xlO-^ 


5.8 xlO-^ 



Note. — () (cr) represents a time average (temporal standard deviation) in the turbulent state. The time average is taken over the quasi-steady 
turbulent state which lasts from fa 3 — 10 orbits at ro depending on the simulation. All RM results are also averaged over a lOrg beam. 

^6 - 90° 

^ 30° 



^The fiducial run. 



Table 2 

RM FOR Initial Field Configuration B 





Resolution 




(Coh> » 


cr(Coh) 






(Coh) 


CT(Coh) ^ 




O-(flMcode) 


200 rg 


120 X 88 


0.0016 


-0.17 


0.27 


-5.9 xlO"* 


9.1 xlO"*' 


0.45 


0.29 


7.7 xlO"^ 


9 xlQ-^ 


100 rg 


120 X 88 


0.0031 


-0.1 


0.25 


-3.4 xlO^* 


8.8 xlO"" 


0.64 


0.25 


6.9 xlO-^ 


4.6 Xl0~^ 


100 rg 


240 X 176 


0.0061 


-0.2 


0.21 


-3.4 xlO"* 


3.8 xlO"" 


0.15 


0.45 


3.2 xlO~^ 


lo-* 


40 r. 


120 X 88 


0.0059 


-0.29 


0.17 


-0.0012 


7.2 xlO"'' 


0.47 


0.33 


1.2 xlO"-'' 


1.9 xlO~^ 


40 rg 


240 X 176 


0.0095 


-0.01 


0.12 


-7.1 xlO"" 


4.6 xlO"'' 


0.13 


0.35 


10"^ 


2.6 xlO-^ 


NOTE.- 


— () ((t) represent 


s a time 


average (tern 


iporal standard deviation) in 


the turbulent 


state. The 


time average 


is taken over 


the quasi-steady 



turbulent state which lasts from 3 — 10 orbits at ro depending on the simulation. All RM results are also averaged over a lOr^ beam. 
= 90° 
'■e = 30° 
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Fig. 1. — Magnetic field unit vectors for initial 
field configuration A. The initial magnetic field is 
confined within the dense torus and is along con- 
stant density contours. 
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Fig. 2. — The mass accretion rate through the in- 
ner boundary {r^i^ = "^rg) for the fiducial run. 
There is large variability in the accretion rate 
on timescales comparable to the orbital period at 
small radii. 
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Fig. 3. — Contour plots of density (p, left) and 
(3 = Sttp/B"^ (right) in the quasi-steady turbu- 
lent state for the fiducial run (averaged from 6 to 
20 orbits at lOOr^). Arrows indicate the average 
magnetic field direction. Dense high (3 plasma is 
confined to a relatively thin disk; the density de- 
creases significantly towards the poles, where the 
plasma is strongly magnetized (low (3). The in- 
ner (rin ~ lO^g) and outer (rout = 120rg) radii 
used in our rotation measure integrals are shown. 
Also shown is a fiducial "beam" with a diameter 
of i'bcam — lO^g a.t 9 = 30° used to calculate the 
rotation measure for a finite sized source. 
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Fig. 4. — Time averaged RM (solid line) and stan- 
dard deviation in time (dashed line) along indi- 
vidual lines of sight (9 = constant) in the turbu- 
lent steady state for the fiducial run (ri„ = lOr^, 
''out — 120rg). Temporal fluctuations in RM in 
the turbulent disk (~ 10° about the midplane) 
are much larger than the average RM. Only in the 
polar regions (~ 40° about the poles) are the fluc- 
tuations in RM smaller than the mean. 
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Fig. 5. — Time averaged RM (solid line) and stan- 
dard deviation in time (dashed line) in a lOrg di- 
ameter beam as a function of viewing angle for the 
fiducial run (rin = lOrg, rout = 120rg). Fluctua- 
tions near the midplane of the disk [9 ~ 90°) are 
somewhat smaller than for the individual lines of 
sight shown in Figure ID Nonetheless, only near 
the poles (9 ~ 30°) are the fluctuations smaller 
than the mean. 




Fig. 6. — dRAIcodc/dlnr for a lOrg beam for 
the fiducial simulation. Time average (solid line) 
and standard deviation in time (dot-dashed line) 
for an equatorial viewing angle and time aver- 
age (dashed line) and standard deviation in time 
(dotted line) for a 30° viewing angle are shown. 
Also shown is an r~^/^ fit to (dRMcodc/dlnr) for 
9 = 30°. The curves are oscillatory in part be- 
cause dRAIcodc/dlur requires taking a numerical 
derivative. 




Fig. 7. — Beam-averaged (fbcam — lOrg) RM 
(solid line) and standard deviation across the 
beam a±{RM) (dashed line) as a function of time 
for the fiducial simulation; one orbital period at 
100 rg is w 3(Mbh/4 x 10^ Mq) days. Note that 
a±{RM) in this Figure is a standard deviation ac- 
cross the beam, rather than a standard deviation 
in time, which is shown in the remaining Figures. 
Top: 9 = 30°; Bottom: 9 = 90°. The emission is 
beam depolarized when (t±{RM)X^ > 1. Because 
the RM at polar viewing angles is primarily pro- 
duced at small radii (Fig. [6|), the characteristic 
timescale for RM and cr±{RM) variability scales 
roughly as rf^^; these calculations have r^n = lO^g. 
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Fig. 8. — Initial magnetic field unit vectors for 
field configuration B. In contrast to the initial 
magnetic field for the fiducial simulations (Fig. [1]) , 
in this case the field is symmetric about the equa- 
torial plane and there is a net inward radial field 
in the equator. 
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Fig. 10. — Density in the equatorial plane as a 
function or radius, averaged over the turbulent 
state for initial field configuration A (solid; Fig. 
[1]) and B (dashed; Fig. [8]). The density is aver- 
aged over a 12° wedge centered on the midplane. 
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Fig. 9. — Contour plots of density (p, left) and 
P = Sirp/B^ (right) in the quasi-steady turbu- 
lent state for initial field configuration B with 
ro — lOOrg and resolution 240 x 176 (averaged 
from 1.5 to 3.5 orbits at lOOr^). Arrows indicate 
the average magnetic field direction. The disk is 
less dense, more strongly magnetized (/3 ~ 1), and 
much thicker than for the fiducial simulation (see 
Fig. [3|). The magnetic field in the corona is less 
ordered than for the fiducial simulation, though 
there is still very low (3 plasma at the poles. The 
inner (nn = lO^g) and outer (rout = SOr^) radii 
used in our rotation measure integrals are shown. 
Also shown is a fiducial "beam" with a diameter 
of Dboam = lOr-g at 6* = 30° used to calculate the 
rotation measure for a finite sized source. 
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Fig. 11. — Plasma (3 in the equatorial plane as 
a function or radius, averaged over the turbulent 
state for initial field configuration A (solid; Fig. [1]) 
and B (dashed; Fig. [8]); /3 is averaged over a 12° 
wedge centered on the midplane. The plasma f3 
is significantly smaller for configuration B, likely 
because of flux freezing of the initial radial field. 
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Fig. 12. — Time averaged RM (solid line) and 
standard deviation in time (dashed line) in a 
Wrg diameter beam as a function of viewing an- 
gle for initial field configuration B (for resolu- 
tion 240 X 176, Tin = lO'^gi ''out = 60rg, and 
ro = lOOrg). The fluctuations in RM exceed the 
mean at essentially all angles. 
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Fig. 13. — dRMcodc/d\nr for a lOr^ beam for ini- 
tial field configuration B (for resolution 240 x 176, 
Tin = lOrg, Tout = 60rg, and ro = lOOr^). Time 
averaged (solid line) and standard deviation in 
time (dot-dashed line) for an equatorial viewing 
angle and time average (dashed line) and standard 
deviation in time (dotted line) for a 30° viewing 
angle are shown. dRMcodc/dlnr decreases with 
increasing radius for the 30° viewing angle while it 
increases with increasing radius for the equatorial 
viewing angle. The curves are oscillatory in part 
because dMcodc/d\nr requires taking a numerical 
derivative. 
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Fig. 14. — The effect of magnetic field reversals on 
the rotation measure is quantified by the Coher- 
ence (see eq. [5] ) . Here Coh is shown as a function 
of time for three different resolutions for the fidu- 
cial run: 120 x 88 (solid hue), 60 x 44 (dashed 
line), and 240 x 176 (dot-dashed line). An equa- 
torial viewing angle and a lOr^ beam is assumed 
for these calculations. 
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Fig. 15. — Coherence (see eq. [3]) as a func- 
tion of time for initial field configuration B at two 
resolutions (with r-m — lOr^, rout = 60rg, and 
ro = lOOrg): 120 x 88 (solid line) and 240 x 176 
(dashed line). An equatorial viewing angle and a 
lOrg beam is assumed for these calculations. Note 
that the Coh is positive at some times, demon- 
strating that RM in the midplane can change sign 
even though there is an initially radial field thread- 
ing the equator of the torus. 



19 



